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Abstract 

In apple {Malusxdomestica Borkh), as in many fruiting crops, fruit maintenance vs abscission is a major criteria for 
production profitability. Growers routinely make use of chemical thinning agents to control total fruit load. However, serious 
threats for the environment lead to the demand for new apple cultivars with self-thinning properties. In this project, we 
studied the genetic determinism of this trait using a progeny derived from the cross between the hybrid INRA X3263, 
assumed to possess the self-thinning trait, and the cultivar 'Belrene'. Both counting and percentage variables were 
considered to capture the fruiting behaviour on different shoot types and over three consecutive years. Besides low to 
moderate but significant genetic effects, mixed models showed considerable effects of the year and the shoot type, as well 
as an interaction effect. Year effect resulted mainly from biennial fruiting. Eight Quantitative Trait Locus (QTL) were detected 
on several linkage groups (LG), either independent or specific of the year of observation or the shoot type. The QTL with 
highest LQD value was located on the top third of LGIO. The screening of three QTL zones for candidate genes revealed a 
list of transcription factors and genes involved in fruit nutrition, xylem differentiation, plant responses to starvation and 
organ abscission that open new avenues for further molecular investigations. The detailed phenotyping performed revealed 
the dependency between the self-thinning trait and the fruiting status of the trees. Despite a moderate genetic control of 
the self-thinning trait, QTL and candidate genes were identified which will need further analyses involving other progenies 
and molecular investigations. 
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introduction 

Organ abscission is a natural process tiiat allows plants to 
remove damaged, senescent or mature organs. It results from the 
development of abscission zones in each organ, even though only 
one zone is activated at each specific developmental stage [1], [2], 
[3] . Fruit abscission has been particularly studied in a number of 
species such as tomato, grape, stone and pome fruits, because of its 
importance in determining fruit crop quantity and quality. In 
apple, fruit abscission occurs at three particular developmental 
stages, first a few days after anthesis, second in June before the 
beginning of exponential fruit growth, and third before ripening 
[4], [3]. Because flowers and fruits are formed in clusters located 
on terminal positions of the shoots, fruit drop involves, in addition 
to competition among inflorescences, and between inflorescences 
and vegetative shoot growth, a competition among developing 
fruits within a cluster [5], [6]. This competition has been described 
as a consequence of the relative position of the fruits within the 
cluster, with the terminal flower (also called "king flower") being 
dominant [7] [3] [8]. The nutritional status of the young fruits, 
through the level of sucrose in the pedicel [9] , as well as auxin and 



GA regulation and transport [3] [8] have been considered as 
factors involved in young fruit drop. In recent experiments, the 
molecular signatures related to fruit abscission induced by thinning 
chemicals have confirmed the involvement of a cross-talk between 
the nutritional status of the fruit and hormonal signalling in 
abscission zone activation [8] [10] [11] [12]. According to these 
authors, unfavourable nutritional conditions and sugar availability 
perceived by the young developing fruits induce simultaneously an 
up-regulation of ABA and ethylene with a down-regulation of GA 
signalling pathways. Thus, we hypothesize that within inflores- 
cences lateral fruitiets may develop poorly due to unfavourable 
nutritional conditions, and because of the reduced sink they 
represent, their hormonal production and perception may be 
altered. This change in hormonal balance may in turn activate the 
development of an abscission zone, as well as a number of cell-wall 
degradation enzymes, such as cellulase, polygalacturonase or 
glycolases [11] [13] [14] [15], leading to fruitlet drop several days 
later. 

In fruit tree industry, considering the huge amount of flowers or 
inflorescences that a fruit tree can bear, fruit load control has 
received particular attention. Indeed, an excess of fruits with 
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respect to vegetative growth may lead to low fruit size and to 
irregular or biennal bearing in many perennial crops, particularly 
in apple, pear, plum, olive, and Citrus [16]. Thinning methods are 
thus widely used to promote fruit abscission and control fruit load 
[17], [18], [19] [1]. In apple, chemical thinning is commonly 
applied up to 30 days after full bloom, this period being considered 
optimal because fruit to fruit trophic competition and the 
detrimental effect of fruit on floral initiation are still low. However, 
the effect of thinning treatments is uncertain and largely depends 
on the cultivar and environmental conditions. Moreover, thinning 
agents such as the benzyladenine (BA) or the Naphtaleneacetic 
acid (NAA) may present a threat for the environment and their use 
is being restricted. This leads to the demand for alternative 
strategies among which the selection of new cultivars with self- 
thinning properties. 

From a genetic point of view, most of the processes involved in 
yield determination vary greatiy among cultivars. Varietal 
differences have been reported in the floribundity and propensity 
to regular bearing in the apple tree [20] [21] [22]. Spurs mutants 
with short internodes and lateral branching have an enhanced 
precocity and a greater fruit-setting ability [23], [24], even though 
they tend to be alternate bearing [25] . Additional evidence of the 
impact of vegetative shoot length on terminal bud floral induction 
and fruit setting has been provided, supporting the assumption of 
within tree variation of fruit set [5] [26]. Describing the genetic 
variability of architectural traits among cultivars and progenies, 
[20] noticed several cultivars and hybrids that exhibited a fruit self- 
thinning behaviour. These genotypes were mentioned as naturally 
maintaining one fruit per inflorescence after fruit set and were 
introduced as parents in several crosses at both Bordeaux and 
Angers INRA stations (France). Despite this strong hypothesis of 
genetic transmission of the self-thinning trait, no study was 
performed so far to confirm its genetic determinism. Natural 
variation of fruit abscission has been reported [27], but this study 
focused on the last developmental stages before ripening rather 
than on fruit abscission just after anthesis. In the present study, in 
addition to genetic control, we hypothesized a possible influence of 
crop load on the ability of the tree to self-thin the fruits because of 
the role of nutritional status of fruits in the natural drop. We thus 
investigated the genetic determinism of the self-thinning trait with 
respect to (i) its stability over years and annual variation of tree 
fruit load and (ii) within-tree trait variation resulting from different 
axis types bearing fruits. A detailed phenotyping of a progeny 
derived from an INRA hybrid and a QTL analysis were 



performed in order to identify the genomic regions involved in 
the genetic variation of this trait. 

Materials and Methods 

Plant Material 

The Fl progeny under study is derived from a cross between the 
hybrid X3263 and the cultivar 'Belrene'. The parents were initially 
chosen for their contrasted architectural traits, 'Belrene' exhibiting 
an erected tree habit (type II according to Lespinasse's classifica- 
tion [28]) while X3263 is considered to have an intermediate 
growth habit (type III). This hybrid was bread at the INRA station 
of Angers and derived from a cross between 'Red Winter' and 
X3177, the latter being itself a hybrid derived from a cross 
between 'Idared' and 'Prima'. X3263 hybrid was described as not 
sensitive to alternate bearing and exhibiting self-thinning trait (Y. 
Lespinasse, personal communication). 

The segregating population is composed of 324 trees, of which 
50 were randomly selected to produce rephcates. Trees were 
grafted onto 'Pajam 1' apple rootstocks and planted in 2005 at a 
single location, the MelgueH INRA Montpellier experimental 
station. All hybrid trees as well as the parents and grand-parents 
present in the orchard were phenotyped, and 271 genotypes, 
including those with two replicates, were used for the linkage map 
construction. Because of alternate bearing behaviour, the number 
of trees observed varied depending on the year. Thus, 286 hybrid 
trees were observed in 2008, which corresponded to the first year 
of flowering of the progeny, 296 were observed in 2009 and 285 in 
2010. 

Phenotyping 

Fruit set was recorded on inflorescences born on two branches 
per tree. These branches, located along the trunks, had developed 
in the same year and were chosen as comparable as possible in 
terms of development. Along each branch, the successive years of 
growth were identified and the shoots born along those sections 
were classified in three types, depending on the length of growth 
units (GUs) that composed them as previously described in [29]: 
shoots were considered short when the length of all their GUs was 
less than 5 cm; they were classified as medium when at least one 
GU was more than 5 cm but less than 20 cm long; and shoots 
were considered long when at least one GU was more than 20 cm 
long. 




Figure 1 . Mean number of inflorescences per shoot for eacKi slioot type and year of observation in a F1 apple progeny derived from 
the cross 'X3263'x'Belrene'. Ax: Axillary Inflorescences; S: Short shoots; M: Medium shoots; L: Long shoots. 
dol:1 0.1 371/journal.pone.0091 01 6.g001 
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Figure 2. Distribution of the mean number of fruits per inflorescence with fruit set depending on the year of observation in a F1 
hybrid apple tree population derived from the cross X3263x'Belrene'. (See also the mean values observed each year for the two parents 
and grand-parents in Table 1). 
doi:10.1371/journal.pone.0091016.g002 



Mid-April of the first observation year (2008), near full bloom, 
the number of inflorescences was counted along each selected 
branch, noting the bearing shoot category. The inflorescences that 
were axillary along the one-year-old shoot of the branch were also 
counted and classified as 'axUlary'. In June, about 70 days after fuU 
bloom (dafb) and after fruit set and natural thinning, the number 
of fruits per inflorescence that had set was counted. The following 
year (2009), the total number of inflorescences per shoot was 
recorded for the second time, at the same period as in the first 
year. The number of fruits per terminal inflorescence was counted 
twice: a first counting was performed in May (about 30 dafb) and 
the second in June (about 70 dafb). In the third year (2010), after 
having checked that the total number of inflorescences observed at 
30 dafb and 70 dafb were highly correlated, the number of fruits 
set per inflorescence was recorded in May only (30 dafb). 

Finally, we obtained a dataset that contained the following 
variables for three years: the total number of inflorescences per 
shoot type along a branch (NIn), the fruit set per branch that was 
estimated as the ratio between these two variables (NIn_s/NIn), 
and the number of fruits per terminal inflorescence (NFI). An 
additional variable was considered for each shoot type: the 
percentage of inflorescences with 1 fruit (%In_lfr) which 
corresponds to its ratio with respect to the total number of 
inflorescences (NIn_lfr/NIn). 



Data Analyses 

All statistical analyses were performed using R software v3.0.2 
[30], with lme4 package for mixed linear model estimation. 
Asreml software was used for the calculation of confidence interval 
on heritabUity values [31]. 

The normality of each variable distribution was checked. When 
not verified, the corresponding variable was transformed with a 
square root transformation before modelling (these variables are 
distinguished in Tables by the star symbol). 

On all variables observed over three years, significance of the 
year, shoot type and genotype effects, and their first order 
interactions were estimated using the following model: 

Pijk = + Gi + Yj + Sk + (G, X Yj) + (G/ X 5/,) + £,;/t, 

where P is the phenotypic value of the shoot k of the genotype i in 
year J, |i is the total average of the population, Gi is the effect of the 
genotype i, Yj is the effect of the year J, Sj^ is the effect of the shoot 
type k, Gi X Yj is the interaction between the genotype and year, 
G, X Sk is the interaction between the genotype and shoot type, 
and Syi; is the residual error. 

Mixed linear models were built for each variable, considering 
the year (Y) and shoot type (S) as fixed effects, the genotype (G) 
and the interactions with the genotype as random effects. The 
models were estimated with restricted maximum likelihood 
(REML) estimation method and effects to be included were 
selected on the basis of AIC and BIG (Akaike and Bayesian 



Table 1. Mean number of fruits per inflorescence with fruit set for each parent (Belrene and X3263) and two grand-parents (Red 
Winter and X3177), depending on the year. 






2008 


2009 


2010 


All years 


Belrene 


2.92 


3.49 


2.17 


2.87 


X3263 


1.82 


2.01 


1.53 


1.77 


Red Winter 


1.80 


2.47 


2.76 


2.31 


X3177 


1.03 


1.36 


1.32 


1.25 



doi:1 0.1 371 /journal.pone.0091 01 6.t001 
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Table 2. Spearman phenotypic (below diagonal) and genotypic (above diagonal) correlation coefficients between the total 
number of inflorescences per branch (NIn), the fruit set, the number of fruit per inflorescence in terminal position (NFI), the 
percentage of Inflorescences with 1 fruit (%ln_1fr). 







NIn* 


fruit_set 


NFI* 


%ln_1fr 


NIn* 


1 


0.02 


0.08 


-0.12 


fruit_set 


-0.47 


1 


0.53 


-0.52 


NFI* 


0.02 


0.37 


1 


-0.87 


%ln_lfr 


0.11 


-0.35 


-0.86 


1 



Variables transformed with a square root transformation are indicated by * symbol. Significant coefficients are indicated in bold. 
doi:10.1371/journal.pone.0091016.t002 



Information Criterion, respectively) minimization. For each trait, 
when G and interaction effects were included, the heritabUity 
of the mean genotype (h^) was estimated as the ratio between 
the genotypic variance and the total variance. Variables were 
considered as heritable if their h^ value was greater than 0.2 
[32]. When an interaction G, x Yj was included in the 
model the heritabUity was calculated as follows: = 

,where n is the number of repetitions per 



genotype (2, in the present case), and n^ is the number of years (3 
in the present case). 

When the interaction G, x Sk was also included, the calculation 
was: 



(4 + 



+ 



2 

GxS 



where n is the number of repetitions per genotype, n,^ is the 
number of years and n^ is the number of shoot types (4 in the 
present case, corresponding to long, medium, short axillary shoots 
along the branches and inflorescences directly inserted in axillary 
positions along the 1 -year-old section of the same branches). 

Confidence interval for heritabUity values were estimated with 
asreml and deltamethod procedure. 

QTL Detection 

Both parental genetic linkage maps and an integrated map of 
X3263 x'Belrene', previously developed [33], were used for QTL 
detection. These maps were constructed using 271 individuals with 
83 SSR markers [34], [35] and 128 SNP markers [36]. A total of 
211 genetic markers were mapped on the X3263 x'Belrene' 
integrated genetic map, which covers 1068 cM over 17 LGs [33]. 
However, markers with similar alleles between parents (hkxhk) 
were discarded for QTL detection, and the map considered 
includes 186 markers. 

QTL analysis was performed using MapQTL6 [37] on the 
mean genotypic values for all variables and on BLUP extracted 
from the linear mixed model for normally distributed variables 
only. QTLs were detected using the interval mapping (IM) (step 
size 1 cM) and multiple QTL mapping (MQM) functions. A QTL 
was declared significant if the maximum LOD score exceeded the 
genome wide LOD threshold calculated over 1000 permutations, 
with a mean error rate of 0.05. When a QTL was declared 
significant after IM, the nearest marker to the LOD peak was 
selected as co-factor for MQM. Each QTL was characterized by 
its LOD score, the percentage of phenotypic variation explained, 
and its confidence interval that corresponded to a LOD score drop 



of 1 or 2 on either side of the likelihood peak. The allelic effects 
were estimated for female and male additivity, and for dominance 
(OUivier, 2002). Female {AS) and male (Am) additive effects 

were computed as Aj = and Am = 



respectively where n^^, fi^ 



are estimated phenotypic means associated with each of the four 
possible genotypic classes deriving from a <abxcd> cross. 
Additive effect represents the contrast between mean value per 
class including a vs b allele inherited from the female parent or c vs 
d allele inherited from the male parent, respectively. Dominance 

eliect was computed as 1)= and corre- 

sponds to an interaction effect between female and male alleles. 
Length of SSR alleles (in bp) at loci associated with putative QTLs 
is indicated in Table SI. 

QTLs were graphically displayed as bars next to the LG on 
which they were identified using MapChart version 2.0 [38]. 
When several QTLs were detected for a trait, a global model 
including all cofactors and their order 2 interactions, considered as 
fixed effects, was built to test for epistatic effects. The significant 
effects were selected on the basis of AIC and BIG criteria 
minimization. This modelling step allowed to estimate the global 
percentage of phenotypic variation (global R^) explained joindy by 
aU the QTLs. 

Search for in silico Candidate Genes and miRNA 
Positioned under QTL Regions 

Among the QTLs detected, we selected those which had a 
confidence interval less than 5 cM for a more in depth 
investigation of putative candidate genes. For this, markers 
flanking the QTL region were first retrieved from the reference 
apple genome, available at Genome Database for Rosaceaae 
(GDR) website (http://www.rosaceae.org/cgi-biii/gdr, 2013 Nov 
1), and their position downloaded from Malus xdomestica whole 
genome vl.O Assembly & Annotation fde "Malus_x_domestica. 
V 1.0. markers.xls". The sequences of markers not listed in this file 
were blasted against Apple genome vl.O contigs with the GDR 
tool NCBI BLAST to recover their positions. The list of predicted 
apple genes located between the positions of the markers flanking 
the QTLs was extracted. 

From the list obtained for each QTL region, we focussed 
particularly on transcription factors and genes which putative 
function suggest a role in development of floral organs, carbohy- 
drates and water supply to the young fruits, and abscission 
processes. We also investigated the possible involvement of 
miRNA genes that could interact with TF or genes of interest in 
the selected QTL zones. For this, the genes sequences were 
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submitted to the 'Preloaded small RNAs/user-submitted tran- 
scripts' tool of psRNATarget server [39] (http:/ /plantgrn. noble. 
org/psRNATarget/?function = 2) against the 207 sequences of 
miRBase (www.miRBase.org, Release 19, August 2012). 

Results 

Trait Variation Depending on the Observation Year and 
Shoot Type 

Among the three years of observation, 2008 corresponded to the 
first year of fruit production. In the following year (2009), most 
trees exliibited a heavy crop load whereas the crop load was lower 
in 2010 (Figure 1). This behaviour corresponds to an irregular 
bearing, characterized by an "on" year in 2009 and "off' years in 
2008 and 2010. The total number of inflorescences observed on all 
the sampled branches varied considerably between years, with 
about 10.000 inflorescences in 2008, 34.000 in 2009 and 6.000 
inflorescences observed only in 2010. The alternate bearing was 
also apparent at the local scale, i.e. the mean number of 
inflorescences per shoot type was higher in the "on" year than 
in the "off' years on three shoot types (Figure 1). The difference 
between years was more pronounced on long shoots than on 
medium and short shoots, this latter being proportionally less 
represented in the "on" year (data not shown). However, the 
difierence in mean number of inflorescences per shoot type 
between "on" and "off' years was not observed for inflorescences 
born in axillary position. In that case, the mean number of 
inflorescences slightly decreased over the three years, indepen- 
dently from the "on" or "off' status of the year. 

Considering the inflorescences with fruit set only, the distribu- 
tion of the mean number of fruits per inflorescence was L-shaped 
for the two "off" years, with a majority of inflorescences 
supporting one fruit and few inflorescences with a high number 
of fruits (Figure 2). However, we found more inflorescences with 
two fruits than a single fruit in the "on" year. As expected, 
'Belrene' had a higher mean number of fruits per inflorescence 
than X3263, whatever the year (Table 1, Figure 2). It is 
remarkable that the parent X3263 had slightly more fruits per 
inflorescence than its own parent X3177 which number of fruits 
per inflorescence was close to 1 whatever the year and its bearing 
status ("on" or "off'). For all parents, grand-parents and 
progenies, the mean number of fruits per inflorescence was higher 
in the "on" year (2009) than in the two "off' years (2008 and 
2010). 

Spearman correlation coefficient values indicated that the 
studied variables were correlated to each other (Table 2). The 
total number of inflorescences (NIn) was poorly correlated with the 
other variables. Interestingly, the slightly negative phenotypic 
correlation with fruitset (—0.47) was almost null for genotypic 
correlation, suggesting that antagonism between number of 
inflorescences and fruitset was entirely environmentally driven. 
Correlation between fruitset and the number of fruits per 
inflorescence was positive (0.37 and 0.53 for phenotypic and 
genotypic correlations respectively) whereas negative values were 
observed between fruitset and the percentage of inflorescences 
with one fruit (—0.35 and —0.52 for phenotypic and genotypic 
correlations respectively). The increase in absolute value when 
considering genotypic correlations is noticeable. Finally, the 
number of fruits per inflorescence was highly and negatively 
correlated to the percentage of inflorescences with 1 fruit (—0.86 
and 0.87 for phenotypic and genotypic correlations respectively). 
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Table 4. Estimates of year (Y), shoot type (ST) and their first-order interactions in the models selected for the total number of 
inflorescences per branch (Nln), the fruit set, the number of fruits per inflorescences (NFI), the percentage of inflorescences with 
one fruit (%ln_lfr), in a F1 apple progeny derived from X3263x'Belrene' cross. 



Year _ shoot type 


Nln* 


Fruit set 


NFI* 


%ln_1fr 


2010 _ Medium 


8.11 


215.92 


75.90 


51.98 


2008 


-2.55 


-31.04 


-20.16 


14.24 


2009 


5.65 


23.66 


25.18 


-8.40 


Axillary 


-3.34 


-38.15 


-23.14 


18.86 


Short 


0.43 


-10.89 


-5.35 


6.78 


Long 


7.11 


22.38 


17.38 


-17.72 


2008 _ Axillary 


5.33 


-24.20 


3.39 


-1.46 


2009 _ Axillary 


-6.72 


26.95 


-4.18 


-3.59 


2008 _ Short 


2.93 


-2.38 


1.54 


-3.64 


2009 _ Short 


-4.01 


4.01 


-7.85 


-0.06 


2008 _ Long 


-9.41 


19.43 


3.29 


3.36 


2009 _ Long 


14.05 


-25.85 


9.49 


4.33 


For estimates, medium shoots 


n year 2010 is the reference. Variables transformed with a square 


root transformation are 


indicated by * symbol. 



doi:1 0.1 371 /journal.pone.0091 01 6.t004 



Genetic and Fixed Effects 

The mixed linear model with three factor.s revealed a significant 
random genetic effect for all the studied variables, except the total 
number of inflorescences per branch (Table 3). For all variables, 
both the year and shoot type and their interaction were included 
as frxed factors in the selected models. The relative influence of 
frxed effects can be analysed through the effect estimates on the 
mean values, which are expressed with the medium shoots in 2010 
as reference (Table 4). This analysis confirms that Nln and NFI 
variables had lower mean values in 2008 and 2010 than in 2009, 
whereas the year effect was inverse on the percentage of 
inflorescences with 1 fruit (%In_lfr). Long shoots had higher 
mean values of total number of inflorescences (Nln), number of 
fi'uit per inflorescence (NFI) and fruit set than medium shoots, 
whereas inflorescences in axillary position had lower mean values 
for these variables. Interactions between years and shoot types also 
revealed that the total number of inflorescences (Nln) was lower in 
the "on" year (2009), compared to "off" years (2008 and 2010) for 
inflorescences in axillary position and located on short shoots 
whereas the inverse was observed for long shoots. The impact of 
these interactions on the mean values of the other variables was 
lower, even though significant. 

The selected models also included interactions between the 
genotype and the year and the shoot type factors for all the studied 
variables. Medium heritabUity values were obtained for the 
number of fruits per inflorescence and the percentage of 
inflorescences with one fruit with confidence interval between 
0.48-0.56 and 0.52-0.62, respectively (Table 3). Fruit set exhibited 
a higher heritability with a confidence interval from 0.61 to 0.70. 

QTL Detection 

Two QTLs were identified on the consensus genetic map for the 
mean genotypic values of the total number of inflorescences per 
shoot type or year (Table 5, Figure 3). A QTL located on LG12 
was detected for the total number of inflorescences in 2008 
(NIn_08), consistently with the significant GxY interaction 
(Table 3). This QTL explained 1 1 % of the variance, and was 
characterized by a male additive effect, and was also detected on 
the male 'Belrene' parental map (with a LOD score of 3.4). A 



second QTL was detected for the total number of inflorescences 
per medium shoot (NIn_M) on LGl, consistently with the 
significant GxST interaction (Table 3). It explained 7% of the 
trait variance and was characterized by a female additive effect. 
This QTL was also identified in the female 'X3263' parental map 
(with a LOD score of 4.1). 

Several QTL were detected for the percentage of the total 
number of inflorescences with one fruit. A close-to-significant 
QTL was detected on LG 1 6 for this variable considered whatever 
the year and the shoot type. It explained a relatively low 
percentage of the variance (6.5%) and was characterized by both 
male and dominance effect. The same genomic region was found 
to control the percentage of inflorescences with one fruit in 2009, 
with similar characteristics. Two other QTLs were detected for 
this trait on LG9 and LGIO, and explained 6.9 and 9.1% of the 
variance, respectively. The global model built with the three 
selected co-factors included a significant interaction between 
CH02b07_XB on LGIO and GDI 42 on LG9, and explained 30% 
of the total variance of the percentage of inflorescences with one 
fruit in 2009. Another QTL located on LGIO near the COL 
marker was detected for the percentage of inflorescences with one 
fruit on spurs. This QTL explained 7.9% of the variance and 
involved mainly a male additive effect even though it was not 
detected on the parental map. 

No or few QTLs were detected for the number of fruits per 
inflorescence whatever the shoot type or the year considered. Only 
one QTL was identified for the number of fruits per terminal 
inflorescence in 2009 (NFI_09) on LG 16, at a similar region than 
previously detected for inflorescences. This QTL explained 8% of 
the trait variability and was characterized by a male additive 
effect. It was also detected on the male parental map with a LOD 
score value of 3.4. By contrast, several QTL were detected when 
the fruit numbers were considered as a percentage, i.e. for fruit set 
related variables (Table 5, Fig. 3). A QTL was detected on the 
distal part of LGIO for fruit set whatever the year or the shoot 
type. It explained 12.4% of the variance and was characterized by 
a male additive and a dominance effect. This QTL was in the 
same genomic region as those detected for the percentage of 
inflorescence with one fruit in 2009 and was detected on 'Belrene' 
parental map with a LOD value of 3.0. Other QTLs were detected 
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Figure 3. Genomic positions of the QTLs detected on the 'X3263'x'Belrene' consensus map. QTLs are represented by boxes, in which 
bold lines represent the LOD-1 confidence interval and extended lines represent the LOD-2 confidence interval. Boxes representing QTLs for the 
number of inflorescences are white, those for the number of fruits per inflorescence are black, for the fruitset traits are pale grey, and for the 
percentage of inflorescences with 1 fruit are dark grey. For trait abbreviations, see Table 5. 
doi:10.1371/journal.pone.0091016.g003 



in the same region for fruit set in 2009 and 2010, and for fruit set 
of spurs. These QTL explained 9.6%, 12.4%, and 17.9% of the 
variance, respectively. The QTL for fruit set in 2009 was 



characterized by a dominant effect whereas that for fruit set in 
2010 was characterized by both a male additive and dominance 
effects. The QTL for fruit set of spurs was characterized by a male 
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additive effect and was also present on the parental map (LOD 
4.0). Another QTL was detected for this trait on LG15, and 
explained 8.1% of the variance. Both QTL explained 13% of the 
variance of fruit set on spurs without interaction between the two 
corresponding co-factors (CH02b07_XB and ch02c09). For fruit 
set in 2008 a QTL was identified on LG12.It explained 12.6% of 
the variance and was characterized by both male and female 
additive effects. For fruit set on long shoots a QTL was located on 
LGl 1, and explained 9.1% of the variance. These two QTLs were 
also identified in the female parental map with a LOD score of 3.5 
and 2.6 respectively (the genome wide threshold being 2.5). 

Screening for Putative Candidate Genes 

Three QTL regions were selected and screened for putative 
candidate genes (gene list in Table 6 and Table S2). Two regions 
were located on LGIO and one on LG9. During this investigation, 
particular attention was paid to transcription factors and predicted 
genes involved in vascular tissues development and sugar 
transport, both of which processes may contribute to the 
maintenance or abscission of fruitlets (Table 6 and Supplementary 
Tables). 

The first genomic region investigated on LGIO was comprised 
between die markers CH02b07 and GD_SNP01867. It spanned 
the length of about 2 Mb and included 243 predicted genes. 
Within this zone, 17 transcription factors (TF) were identified, 
among which three were previously found to be involved in 
vascular tissues development and sugar transport in Arahidopsis 
thaliana. In particular, we identified a Vascular Relattxi NAC- 
domain TF homolog to Arahidopsis thaliana VND7 (see Ath and 
MDP numbers in Table 6). Among the genes located within the 
QTL confidence interval we also identified a UDP-glucosyl 
transferase involved in the biosynthesis of poly-saccharides and 
two carbohydrate trans-membrane transporters (SUC 1 and 
SUC2), the second one with two copies. Even though several 
SUC copies were present on almost all LGs, SUC 1 had only three 
copies and SUC 2 exhibited 6 copies along the apple reference 
genome. 

In the same zone, we also found two copies of SCARCROW- 
LIKE 21 protein (SCL21). SCL genes belong to the large GRAS 
family and are involved in a number of developmental processes 
such as GA responses controlling flowering, shoot and root apex 
development or xylem patterning. A third copy of SCL21 was 
present in the second QTL zone on LGIO, close to CO marker. 
These three copies are the only ones present along the apple 
genome for SCL21 gene. 

The second QTL zone investigated on LGIO was comprised 
between the markers GD_SNP00360 and COL, and spanned 
1.2 Mb. It contained 247 predicted gene sequences, including ten 
TF. Among these, we identified another predicted SCL21 TF 
located close to the COL marker. Together with the two SCL21 
previously identified, these three copies are the only ones present 
along the apple genome. In addition this QTL region contained 
one copy of AGAMOUS and two copies of AGAMOUS-LIKE 12 
(AGL12). Among the four copies of AGL12 found in the apple 
reference genome, two were located on LG8 and the two other 
ones were in LGIO, included in the present QTL zone. We also 
identified a TF homologue to SHOOT MERISTEMLESS 
(STM): five copies were found on the apple reference genome, 
two on LGIO, one copy contained in the QTL zone whereas the 
other one was slightly above. 

Within the QTL confidence interval, we also identified 
WRKY65, a member of the large WRKY family. This family is 
involved in the complex regulation of senescence, especially in 
response to biotic and abiotic stress. Only three copies of 



WRKY65 were found in the reference genome on LG5, LGIO 
and LG15. 

Several genes involved in sugar synthesis were also identified 
within this QTL. These include genes homologue to ATTPS 1 and 
ATTPS7, two trehalose synthases, and IRXl (also known as 
CESA8) involved in cellulose synthase. Seven other copies of 
IRXl were present on the apple reference genome, 6 on LG5 and 
one on LGl 6. 

Three Md-miRNAs were listed in the Predicted miRNA/ 
Target Pairs analysis among the genes included in this LG 1 0 QTL 
confidence interval: Md-miR390 (6 isogenes), Md-miR7124 (2 
isogenes) and Md-miR482a-p (list in Table S3 and S4). The first 
miRNA has a cleavage function on MDPOOOO 158644 gene, an 
homolog of AT3G14840.2 annotated as Leucine-rich repeat trans- 
membrane protein kinase. The two others have a translational 
inhibition effect on MDP0000300617 and MDP0000735861, 
homologues of AT4G 18760.1 a receptor like protein 51 (RLP51) 
involved in signal transduction in plasma membrane and 
AT5G36930.2, a disease resistance protein of TIR-NBS-LRR 
class family involved in signal transduction, defense response and 
apoptosis, expressed in different growth stages, respectively. 

The third QTL region investigated was located at the top of 
LG9. This region spanned 4.8Mb and contained more than 1.100 
predicted genes and 21 TF. Among these we identified a TF 
homologous to SEPALLATA 1 (SEPl) and three copies of 
SEPALLATA2 (SEP2), among the four found in apple reference 
genome, and also known as AGAMOUS-LIKE 4 and 2, 
respectively (AGL4 and AGL2). Another member of the large 
AGL family, AGAMOUS-LIKE 8 (AGL8), also known as 
FRUITFULL, was identified in this genomic region. AGL8 is 
involved in a complex regulatory network in which ASSY- 
METRIC LEAVES (ASl), another TF involved in leave and 
gynoecium patterning plays a major role. One TF homologous to 
ASl was also present in this genomic region. 

Within the confidence interval of the LG9 QTL we identified 
seven copies of the Senescence-Associated-Genes SAGlOl, among 
the 18 copies found in the reference genome, and 3 genes 
homologous to GLIP7 which function is similar to SAG genes. In 
addition, several genes potentially involved in the glucose/ cellulose 
biosynthesis pathways were identified: these include 17 glycosyl- 
transferases and two galacturonases genes. Several auxin trans- 
porters of ABCB family (ABCB 19 and PGP17) and expansins were 
also identified. 

Discussion 

The detailed phenotyping performed in the present study 
provides a global view of the relationships between the bearing 
status of the trees, in relation with alternative bearing behaviour, 
fruit set, and the maintenance versus abscission of young fruits 
within inflorescences in a range of genotypes. Our results confirm 
the interdependence among processes occurring at local (inflores- 
cences) and global (whole tree) scales, which are likely to rely on 
the nutritional status of the young developing fruits. In particular, 
we found that the mean number of fruits per inflorescence was 
higher in the 'on' year, when the trees were in a high crop load 
status (Fig. 1, Table 4). This result may appear counterintuitive 
with respect to the usual acceptance that fruit bud density is 
negatively correlated to fruit set and to the number of fruits per 
spur [40] [41]. However, these statements rely on results obtained 
after artificial manipulation of fruit buds or flowers, followed by a 
readjustment of fruit set within the trees. Rather, in the present 
study, the bearing status of the tree was represented by the total 
number of inflorescences per branch, positively correlated to the 
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mean number of inflorescences with one fi'uit, both these variables 
contributing to the 'on' status of the trees. By contrast, the total 
number of inflorescences per branch was not correlated to the 
mean number of fruits per inflorescence and to the percentage of 
inflorescences with one fruit, these two variables being influenced 
by the genotype (Table 2 and 3). Finally, variables defined as 
percentage, such as percentage of inflorescences with one fruit and 
percentage of fruit set were more appropriate for analysing the 
genetic control of tree fruiting behaviour than variables dependent 
on tree sampling. 

Fruit set and the number of fruits per inflorescence were 
observed on diflerent shoot types. The three categories correspond 
to diflerent number of leaves and leaf areas, and consequently to 
diflerent ratio between number of fruits and number of leaves, 
which has been considered as a key factor for fruit set, 
development and final quality [42] [43] . As found for year factor, 
shoot type had a strong effect on the studied variables. Among the 
shoot types, the lowest year effect was observed on short shoots. 
This shows that this shoot type is the most suitable for highUghting 
genotypic diflerences. It demonstrates that short shoots are less 
sensitive to environmental effects, and thus may be the most suited 
for further study of this character on other progenies or elite 
hybrids under selection. 

From a genetic point of view, the total number of inflorescences 
per branch appeared not afiected by the genotype. This probably 
results from the fact that this number depends on the size of the 
sampled branches rather than on the genotype. As a consequence, 
its intcr('st r(;lic's on its capability to reprc'scnt the bearing status of 
the tree which necessitates counting variables [44] [45]. Other 
variables such as fruit set and percentage of inflorescences with one 
fruit were more impacted by the genotype. Because all these 
variables were not normally distributed, it would have been 
justified to use generalized mixed model [46]. However, such 
models are particularly difficult to parametrize, estimate and 
int(-rprc-t and we finally decided to estimate the genetic efiect with 
a standard linear mixed model. As a consequence, we also 
preferred using the mean genotypic values for QTL detection 
rather than the BLUP, even though QTL detection was performed 
with both types of phenotypic variables. The heritability values 
estimated from these models were relatively moderate when 
compared to values previously obtained for resistance [47], fruit 
quality [48] [49] or architectural variables [50] [51] [52] [53]. 
However, the heritability values estimated here are specific to this 
segregating population, and may not reflect the heritability of this 
trait overall. 

Several QTL were detected along the apple chromosomes. 
Consistentiy with heritability values more QTL were detected for 
fruit set and percentage of inflorescences with one fruit than for 
counting variables such as the number of inflorescences or the 
number of fruits per inflorescence. H()^\'c\cr, tlu^ percentage of 
variation explained by these QTL r(;maincd lo^v when compared 
to lu'ritahility \'alues. This sugg(;sts that tiu; existence of "missing 
heritability" probably resulting from the complex nature of the 
self-thinning traits. Often identified in genome wide association 
(GWA) and genomic selection studies [54] [55], this missing 
heritability highlights the limits of marker assisted selection and the 
difliculties to correctly predict the phenotypes in genomic 
selection. 

Moreover, more QTL were identified for the 'on' year (2009) 
than for the two 'off years (2008 and 2010), and for variables 
measured on short shoots than on medium or long shoots. 
Surprisingly, the allelic effects were more often due to male 
('Belrene') alleles rather than to female (X3263) alleles. Since the 
X3263 parent had less fruits per inflorescence than 'Behene', the 



QTL can be interpret as resulting from a "negative" impact of the 
presence of one out of the two 'Belrene' alleles, which promoted 
more fruit set and less inflorescences with one fruit in contrast with 
X3263 alleles. QTLs with female (X3263) additive effect were 
located on LGl and LG9, whereas the QTLs with male additive 
effect were on LGIO (top and bottom of LG), LGl 2 and LGl 6. 
Several QTLs involved a parental and a dominance effect, or only 
a dominance allelic effect. 

The QTL located at the bottom of LGl and controlling the 
number of inflorescences born on medium shoots overlaps with a 
QTL controlling the vegetative budbreak dates in the same 
progeny [33], and co-locates with QTLs identified in another 
progeny, 'Starkrimson' X 'Granny Smith' for annual yields or 
cumulated yields [44] and for the percentage of branching [56]. 
This zone is also very well-known for resistance traits since it 
contains the RviG (Vf) major scab resistance gene and a major scab 
resistance QTL [57] [58]. However, this zone was not further 
investigated for candidate genes because of its large confidence 
interval. 

The QTL identified on the top of LG9 and controlling the 
percentage of inflorescences with one fruits in the 'ON' year co- 
locaUzed with a major QTL controlling vegetative and floral 
budbreak in the same progeny, as well as the green point' variable 
in both 'Starkrimson' x'Grarmy Smith' and 'X3263' x'Belrene' 
progenies [33]. This zone was also identified as controlling 
vegetative budbreak in other progenies [59]. Among few QTL 
studies that have been performed on fruit abscission on other 
crops, most of them dealt with abscission close to harvest, e.g. [60] 
on rice or [61] on melon. [60] noticed that the QTL zone involved 
in rice grain abscission coincides with a QTL previously detected 
for seed dormancy and that could belong to a domestication- 
related block of genes. 

Candidate Gene Identification 

Two major approaches are commonly used to dissect complex 
and quantitative traits, i.e. genome-wide scanning and candidate 
gene approach. In this study, we identified chromosomal regions 
controlling the quantitative traits and performed a candidate gene 
approach using only those portions of the genome located within 
the confidence interval of three major QTLs. However, despite the 
use of both approaches, the practicability of our candidate gene 
approach is limited by its reliance on the first \ crsion of the apple 
genome. Recent findings indicate that up to 15"(} of contigs, and 
therefore of predicted genes, might be misplaced on the 'Golden 
Delicious' genome. Therefore, our identification of putative 
candidate genes is subject to caution, and might not be exhaustive. 
Furthermore, the candidate gene approach is also limited by its 
reliance on the known or presumed biology of the phenotype 
under investigation. In this study, we focussed our attention on 
genes and TF and miRNA which have been shown to be 
potentially involved in sugar synthesis and transport, auxin and 
GA regulation and transport, as well as fruitlet abscission zone 
activation and development. 

The genomic region identified on LG9 contained several TF 
involved in flower and fruit development, and was particularly rich 
in AGAMOUS-like family (AGL4, AGL2, and AGL8) genes. Both 
AGL2 and AGL4 TF have been recently described for their main 
role in apple fruit development and ripening, especially in 
hypanthium tissues [62]. AGL8 has been shown to be involved 
in the coordinated growth of cell tissues and, interestingly, its 
absence was shown to block sUique elongation after fertilization in 
Arabidopsis thaliana [63]. This region also contains putative can- 
didate genes potentially involved in abscission zone development 
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such as glycosyl-traiisferases, galacturonases, expansines and auxin 
transporters of the ABCB family [64] . 

As most members of Senescence- Associated-Gene (SAG) family, 
SAGlOl is involved in plant defenses and responses to stress, and 
more specifically in responses to wounds and pathogen attacks 
[65] . Induced by sugars, it is mainly expressed in leaves but also in 
flowers in Arabiwpsis thaliam (tair). We could hypothesize that in 
fruitlets with poor nutrient and sugar supply, expression of this 
gene might be activated during the abscission process, making this 
gene a potential marker of fruitlet abscission. 

The two genomic region located on LG 1 0 appeared particularly 
important for fruit set, and for the percentage of inflorescences 
with fruit. The genomic region located at the top of LG 1 0, near 
CH04c06 and CH02b07, was the most important for fruit set 
variables, which all co-localised in this region. This region also co- 
localized with a QTL previously detected for the 'green point' 
variable in a 'Starkrimson' x'Granny Smith' population [33]. 
Within this region we identified a putative candidate gene 
homologous to the TF VND7. In poplar, this TF was shown to 
be implicated in secondary cell wall biosynthesis and to be a 
master regulator of xylem vessel element dilferentiation in poplar 
[66] [67]. We could hypothesize that, in apple, the expression of 
this TF might vary among fruitiets of the same inflorescence, and 
might lead to differences in water supply among those fruitiets. We 
also identified genes homologous to SUCl and SUC2, both 
involved in sugar loading and retrieval in and from the phloem 
[68] . Variation in the expression of these genes among fruitiets of 
the same inflorescence could aflect sugar supply and this, in turn, 
possibly leading to differences in fruitiet growth kinetics and to the 
development of dominant fruits. 

Two genes homologous to SCL21 were also identified. These 
genes belong to the large GRAS family and are known to be 
involved in a number of developmental processes such as GA 
responses controlling flowering, shoot and root apex development 
or xylem patterning [69]. As for VND7, up-regulation of this gene 
could increase the number of conducting vessels, and lead to 
varying levels of water supply. SCL21 has been recentiy described 
in barley as regulated by miRNA171 and interacting with 
Phytochrome PATl [70]. 

The second QTL zone located on LG 1 0 in the region of the 
COL marker is an important genomic region which has previously 
been detected in many studies involving the columnar trait [51] 
[71] [72] [73]. Other studies, which did not involve the columnar 
mutation, also highlighted this zone for other traits such as 
maximum internode length in specific years [54] and time of 
vegetative budbreak [33]. Interestingly, this zone also corresponds 
to a number of QTLs controlling fruiting behaviour and regularity 
in a 'Starkrimson' x'Granny Smith' progeny, including precocity, 
number of seeds per fruit, and Bieniial Bearing Index (BBI) [44] . 
Several TF were identified within this QTL region, such as 
AGL12. In Arabidopsis thaliana, AGL12 was found to contribute 
to tissues and organ diflerentiation, especially floral organs [74], 
and more specifically in late flowering [75]. Recentiy, AGL12 was 
found to be specifically expressed in the abscission zone (AZ) in 
conjunction with SHOOT MERISTEMLESS (STM), also 
present in the zone, which was down-regulated during abscission 
of tomato fi"uit [76]. Both genes might have a strong effect on 
apple fruitiet self- thinning. Among the WRKY family [77], 
WRKY65 has been shown to be involved in senescence triggering 
after sugar starvation [78]. We could hypothesize that this gene 
might be expressed in some fruitiet pedicels following a decreased 
sugar supply, thus triggering senescence of tissues and leading to 
fruitiet drop. Other genes were also identified in this QTL region, 
including two putative trehalose synthetases. Trehalose is a 
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disaccharide formed by a 1,1-glucoside bond between two a- 
glucose units, and is an effective signal molecule that may have an 
essential role in developmental processes [79]. Variation in 
trehalose synthase expression in fruitlets of a same inflorescence 
might thus be important for fruit development. Finally, a gene 
with sequence similarity to IRXl was identified within the 
confidence interval of this QTL. IRXl was shown to be involved 
in secondary cell wall biosynthesis. In Arabidopsis thaliana, the loss of 
IRX 1 function provokes mutants with abnormal xylem formation, 
reduced cellulose content, and enhanced drought and osmotic 
stress tolerance [80]. 

Conclusion 

In this study, we performed a QTL analysis to study the genetic 
determinism of complex traits imolved in the self-thinning 
character of young apple fruits. This genetic dissection allowed 
us to identify moderate to strong QTLs using a Fl segregating 
population, and to screen for putative candidate genes located 
within the confidence interval of three major QTLs. This 
investigation allowed us to identify several genes potentially 
involved in the formation of conducting vessels, the synthesis 
and transport of sugars and hormones, and finally the activation 
and development of the abscission zone. Despite its major 
importance, the self-thinning character of fruitlets is not a 
common trait in apple, and very few elite cultivars have the 
ability to induce fruitlet abscission. Thus, the vahdation of QTLs 
in additional cycles of alternate bearing and in different genetic 
background, as well as the identification of new QTLs controlling 
self-thinning will be of major importance in future studies. Finally, 
we believe that the identification of QTLs and genes responsible 
for variations of this character will have a strong impact in future 
apple breeding programme, and wiU ultimately contribute to a 
major reduction in the use of thinning reagents by growers. 
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